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The lowest eigenmodes of the Dirac operator are related to the dynamical breaking of the chi- 
ral symmetry in Quantum Chromodynamics (QCD). In our work we construct quark propagators 
which exclude a varying number of the lowest Dirac eigenmodes and study the influence thereof 
on meson correlators and the meson spectrum. That procedure partially restores the chiral sym- 
metry (in the valence sector) and we observe degeneracies in the spectrum while confinement 
seems not to be affected. 
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1. Motivation 



The Banks-Casher relation [1] associates the density of the smallest non-zero eigenvalues to 
the chiral condensate. The importance of the lowest Dirac eigenmodes in QCD has been investi- 
gated within several studies: in [2-4] it was shown that the pseudoscalar and axial vector corre- 
lators can be approximated in the medium to large temporal distance region when taking only the 
lowest Dirac eigenmodes as input for the quark propagators into account. In [4 - 7] this was used 
to improve the signal of hadron correlators by averaging the low mode contribution to the quark 
propagator over all possible positions of the sinks on the lattice. The authors of [4, 5] studied the 
influence of low mode removal in the quenched approximation using the overlap operator [8,9]. 

In our study [10] we construct meson correlators out of reduced quark propagators which 
exclude a varying number of the lowest Dirac eigenmodes. Using these propagators we compute 
meson masses and we study the possible recovery of the degeneracies of (would be) chiral partners 
in the spectrum. Moreover this gives us insight on the role of the lowest Dirac eigenmodes for 
confinement and thus on the relation of chiral symmetry and confinement. Our modification is 
applied only to the valence quark sector; the configurations have been obtained with two dynamical 
light quark flavors without truncation. 

The two lightest quark flavors in QCD are much lighter than the typical QCD scale. Neglecting 
the mass of these two light quarks, the Lagrangian is invariant under the symmetry group 



SU(2) L xSU(2) i? xU(l) y xU(l) A . 



(1.1) 



The single flavor symmetry U(l)y conserves the baryon number density. The axial vector part of 
the chiral symmetry SU(2)^ x SU(2)# mixes states with opposite parity; the nondegenerate masses 
of parity partners in nature indicate the dynamical breaking of the axial symmetry whereas the 
isospin symmetry is (approximately) preserved in the vacuum. The single flavor axial symmetry 
U(1)a mixes the currents of the same isospin but opposite parity. The latter is not only broken 
spontaneously but also explicitly by the anomaly. 

We will investigate the influence of low mode 
reduction on the pseudoscalar, vector, scalar and 
axial vector currents with the interpolating fields 
given in Table 1 . In a chirally symmetric world, 
the currents number three and six of the table 
would get mixed via the isospin axial transforma- 
tion and the currents one and five would get mixed 
via the single flavor axial transformation U(1)a- 

We did not include isoscalars in our work 
(these need disconnected terms, notoriously dif- 
ficult to compute), otherwise we could additionally study the possible recovery of degeneracies in 
the spectrum between the h\ and current four, between fo and current one as well as between the 
77 and current five, which all would pairwise mix via the axial vector transformation of the chiral 
symmetry. For a detailed discussion of the mixing of different currents under the chiral symmetry 
see [11]. 
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Table 1: The mesons and the corresponding in- 
terpolating fields studied in this work. 
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2. Method 

Most lattice Dirac operators D (an exception is the overlap operator [8, 9]) are nonnormal 
operators, thus one has to distinguish between left eigenvectors (L,| and right eigenvectors \Rj) to 
a given (complex valued) eigenvalue A,. Then the spectral representation of D reads 

N 

D = £ A,- \Ri) (Li\ . (2.1) 

The hermitian Dirac operator D5 = 75D is normal and has real eigenvalues fa. Alternatively D can 
be written in terms of the spectral representation of D5 in the following way: 

N 

D = £^751^)^1. (2.2) 

i=\ 

We want to calculate the meson spectrum using reduced quark propagators which leave out a 
part of the low lying Dirac spectrum. Therefore we split the quark propagator 5 = D 1 into a low 
mode part (lm) and a reduced part (red), e.g., using the eigenvalues and eigenvectors of D$, 

s = £ Mf 1 h> (vj I 7s + £ AT 1 k> ( v i I 75 = S lm5{k) + 5 red5(fc) . (2.3) 

i<k i>k 

Hence we can obtain the reduced part of the propagator by subtracting the low mode part from the 
full propagator 

Sred5(/t) = S — S lm5 ( k j . (2.4) 

Note that the low mode part has to operate on the same sources as S in order to obtain the correct 
reduced propagators. The corresponding separation using the left and right eigenvectors and the 
eigenvalues of D can be done in an analogous way. 

3. The setup 

For our calculation we used 161 gauge field configurations [12, 13] of lattice size 16 3 x32. The 
lattice spacing is a = 0. 144(1) fm and thus the spatial lattice sizeL = 2.3 fm. The configurations are 
obtained with two light degenerate dynamical quark flavors with an AWI-mass of m = 15.3(3) MeV 
and a corresponding pion mass of m % = 322(5) MeV. For the dynamical quarks of the configura- 
tions as well as for the valence quarks of our study the so-called Chirally Improved Dirac operator 
[14, 15] has been used. This operator is an approximate solution to the Ginsparg-Wilson equa- 
tion and therefore exhibits better chiral properties than the Wilson Dirac operator while being less 
expensive by an order of magnitude in comparison to the chirally exact overlap operator. 

On these configurations we computed the lowest 256 eigenvalues of D and the lowest 512 
eigenvalues of D5, see Fig. 1, using ARPACK which is an implementation of the Arnoldi method 
to calculate a part of the spectrum of arbitrary matrices [16]. 

In Fig. 2 we show histograms of the lowest 256 eigenvalues of D and D5. Additionally in the 
same figure the integrals over these histograms are given. 
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Figure 1: The lowest 256 eigenvalues of D (left) and the lowest 512 eigenvalues of D5 (right) for all 161 
configurations. 
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Figure 2: L.h.s.: histograms of the lowest 256 eigenvalues of D (red) and D5 (blue). R.h.s.: integrated 
histograms. 

4. Results 

Before we study the effect of low mode reduction, we look at meson correlators which are build 
solely out of the low mode part of the propagator, Sww or S im5 ^, respectively (following the ideas 
in [2-4]). In Fig. 3 we show the pion and rho correlators (interpolators one and three of Tablel) 
from the low mode only propagators 5 , lm< - 12 s) an d Sim5(32) in comparison with the correlators from 
full propagators. As can be seen in both cases, four times less eigenvalues of D5 in comparison 
to eigenvalues of D are already sufficient to obtain a similar quality approximation of the full 
correlators. Since the eigenvalues of D5 seem to be more effective for the low mode saturation of 
meson correlators, we will concentrate on the eigenvalues of D5 throughout the rest of this work. 

In Fig. 5 we show the correlators under D 5 low mode reduction in comparison to the correlators 
from full quark propagators. We find that negative parity mesons (currents 1-4) become heavier 
under low mode reduction whereas positive parity mesons (currents 5-6) at first become lighter and 
after subtracting more than 8-16 modes, their mass begins to increase as well. Except for the pion 
from current one, all correlators still show a clear exponential form (at least as good as for the full 
propagators) from which we conclude that we are still confronted with bound particles. The loss 
of the exponential behavior in current one can qualitatively be explained with the strong sensitivity 
of the pion on the lowest Dirac modes in the sense that very few modes are sufficient to saturate 
the pseudoscalar correlator in the medium to large time region. Moreover we find that low mode 
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Figure 3: Low mode contribution to the pion (left) and rho (right) correlators. 



1400 



1200 



1000 



# of subtracted eigenmodes 
1 6 32 64 



128 



r 8oo - 



600 



400 



200 




removal reduces the noise of the correlators which is most obvious for the positive parity currents. 

We performed effective mass fits for all 
the reduced correlators and combined the re- 
sulting mass values, as a function of the re- 
duction level, in Fig. 4. For the pion from 
current one we were only able to perform fits 
for very low reduction levels. The masses 
of current three and six are found to be de- 
generate from truncation level 16 onwards 
which corresponds to having subtracted all 
eigenvalues with a magnitude equal to or less 
than twice the bare quark mass of our set- 
ting. Furthermore current four becomes de- 
generate with currents three and six which 
hints towards a larger degeneracy which also 
might include the h\ which we did not calcu- 
late. Whereas we observe chiral restoration 
in the 7=1 states, we cannot draw any conclusions for J = states since the would-be chiral part- 
ner of the ciq (current five) is the pion (current one) which gets destroyed under low mode reduction. 
Current two which couples to the pion due to PCAC (partially conserved vector current) does not 
mix with the scalar so there is no reason these two states should become degenerate. For further 
discussion see [10]. 
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Figure 4: The change of the masses of the mesons un- 
der low mode removal of the ingoing quark propaga- 
tors. The labeling number corresponds to the currents 
from Table 1 . 



5. Conclusions 



We studied the isovector meson spectrum using valence quark propagators where a varying 
number of the lowest Dirac eigenmodes has been removed. The density of the lowest eigenmodes 
is related to the spontaneous breaking of the chiral symmetry [1]. We find restoration of the chiral 
symmetry in / = 1 states when only 16 eigenmodes have been subtracted whereas the exponential 
form of the correlators does not suffer from this modification. An exception is the pion correlator 
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for which the lowest Dirac eigenmodes are crucial, its exponential behavior rapidly gets lost when 
subtracting low modes. We conclude that the low lying Dirac modes - although providing an 
important signal for spontaneous chiral symmetry breaking - are not important for the binding 
mechanism of quarks at least in the isovector mesons (with the exception of the pion). 




Figure 5: Correlators under low mode reduction. The number in parentheses indicates which current is 
considered, cf. Table 1 . 
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